Stars and Halos of Degenerate Relativistic Heavy-Neutrino and Neutralino Matter 



OO 
On 
On 



Neven Bilic 1,2 , Faustin Munyaneza 1 , and Raoul D. Viollier 1 
1 Department of Physics, University of Cape Town, Rondebosch 7700, South Africa 
2 Rudjer Boskovic Institute, 10000 Zagreb, Croatia 
(February 5, 2008) 

Heavy-neutrino (or neutralino) stars are studied using the general relativistic equations of hy- 
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an interacting astrophysical system are discussed in detail. 
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I. INTRODUCTION 



One of the most tantalizing puzzles of this universe is the issue of dark matter, the presence of which is inferred 
from the observed flat rotation curves in spiral galaxies ^,|^] , the diffuse emission of x-rays in elliptical galaxies and 
clusters of galaxies, as well as from cluster dynamics. Primordial nucleosynthesis entails that most of baryonic matter 
in this universe is nonluminous, and such an amount of dark matter falls suspiciously close to that required by galactic 
rotation curves. However, although a significant component of dark matter in galactic halos is presumably baryonic 
[3 1, the bulk part of dark matter in this universe is believed to be nonbaryonic. Many candidates have been proposed 
^v^j 1 both baryonic as well as nonbaryonic, to explain the dark matter paradigm, but the issue of the nature of dark 

[ matter is still far from being resolved. 

One of the most conservative candidates for nonbaryonic dark matter are, of course, massive neutrinos. In this paper 
we are particularly interested in neutrinos with masses between 10 and 25 keV/c 2 , as these could form supermassive 
degenerate neutrino stars, which may explain, without invoking the black-hole hypothesis, some of the features 
observed around the supermassive compact dark objects with masses ranging from 10 6 5 Mg to 10 9 5 Mq . These 
have been reported to exist at the center of a number of galaxies ]5|-|8| including our own [p|-pd[ and quasistellar 
,^3 ■ objects. It is interesting to note that neutrinos in this mass range can also cluster around ordinary stars, and thus 
Or these neutrinos could account for at least part of galactic dark matter. A further motivation for studying the collapsed 
structures of heavy neutrino matter is the recent increased interest in fermionic cold dark matter models |H| in which 
S-^ | massive neutrinos play an important role in structure formation in the early universe. 

A 10 to 25 keV/c 2 neutrino is in conflict neither with particle and nuclear physics experiments nor with astrophysical 
observations . On the contrary, if the conclusion of the LSND collaboration which claims to have detected P M — > v e 
^ , flavor oscillations ]T7]] is confirmed, and the quadratic see-saw mechanism involving the up, charm, and top quarks 
p8j , p9j , is the correct mechanism for neutrino mass generation, the v T mass may be between 6 and 32 keV/c 2 [ p0| , 
which is well within the cosmologically forbidden range. It is well known that such a quasistable neutrino would 
lead to an early neutrino-matter dominated phase, which may have started as early as a couple of weeks after the 
big bang. Thus, a critical universe that remained neutrino-matter dominated all the time would have reached the 
current microwave background temperature in less than 1 Gyr, i.e., much too early to accommodate the oldest stars 
in globular clusters, cosmochronology, and the Hubble expansion age. 

It is conceivable, however, that, in the presence of such heavy neutrinos, the early universe might have evolved 
quite differently than described in the homogeneous standard model of cosmology. Neutrino stars may have emerged 
in local condensation processes during a gravitational phase transition, shortly after the neutrino-matter dominated 
epoch began. The latent heat produced in such a first-order phase transition, apart from reheating the gaseous phase, 
might have reheated the radiation background as well. Annihilation of heavy neutrinos into light neutrinos via the 
Z a would take place in the interior of neutrino stars |2l],[l6| . Both these processes would decrease the density of heavy 
neutrinos, as perceived today, and also increase the time which photons need to cool down to the present microwave 
background temperature. Thus a quasistable neutrino in the mass range between 10 and 25 keV/c 2 is presumably 
not in contradiction with cosmological and astrophysical observations pq ] . 

In fact, it has recently been shown p2| , p0|p2] | that degenerate neutrino stars Jl6| , pl]j2^j2^ ] may indeed have been 
formed during a gravitational phase transition in the early universe. Whereas the existence of this first-order phase 
transition is firmly established in the framework of the Thomas-Fermi model at finite temperature p^] , the microscopic 
mechanism through which the latent heat is released during the phase transition and dissipated into observable and 
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perhaps unobservable matter or radiation remains to be identified. At this stage, however, it is still not clear whether 
an efficient dissipation mechanism can be found within the minimal extension of the standard model of particle physics 
or whether new physics is required in the right-handed neutrino sector. We therefore have to assume in the following 
that such an efficient dissipation mechanism exists, in order to make sure that fermions can actually settle in the state 
of lowest energy in a time much shorter than the age of the universe. 

In this paper, we focus primarily on gravitationally clustered, degenerate nonbaryonic matter consisting of two 
species of weakly interacting stable or quasistable fermions: one with a mass around 15 keV/c 2 which we subsequently 
call "neutrino" , and the other with a mass around 1 GeV/c 2 which we henceforth call "neutralino" . The chosen 
neutralino mass, although a little bit on the low side, offers the possibility of replacing the neutralino with a neutron, 
as the strong-interaction effects of the neutron in neutron star matter can be simulated by an effective mass. Of 
course, this substitution makes sense only as long as the binding energy of the neutron is larger than the Q value for 
the neutron decay, so that the neutron can be considered stable in neutron star matter. 

It is interesting to note that a variety of similar scenarios can be treated within the same framework: Apart from 
a neutrino halo around a neutron star, one could also study a neutrino halo around a white dwarf or around an 
ordinary star |^3|l , since all these baryonic stars can be approximated using similar polytropic equations of state which 
eventually result in the same nonlinear differential equations of the Lane-Emden type. Moreover, by varying the 
polytropic index of the equation of state, one can also investigate the properties of a cold neutrino star immersed in 
a hot radiation field, or in a hot baryonic background, or in a vacuum with nonzero energy density, which all may 
have played a role in the formation process of primordial neutrino stars. Thus the study of this simple interacting 
neutrino-neutralino system allows us to learn a great deal about the properties of gravitationally clustered baryonic 
and nonbaryonic matter. 

This paper is organized as follows: In Sec. [n] in a general relativistic framework, we discuss the properties and 



implications of degenerate neutrino (and neutralino) stars and their Newtonian and ultrarelativistic limits. In Sec. [II 



we generalize the Tolman-Oppenheimer-Volkoff (TOV) equations to include gravitationally clustered, degenerate 



nonbaryonic matter, consisting of neutrinos and neutralinos. Our results are summarized in Sec. IV 



II. DEGENERATE NEUTRINO STARS 

A spherically symmetric cloud of degenerate neutrino matter can be characterized by its mass density p v {r), pressure 
P v (r), and the metric in the Schwarzschild form p5| 

ds 2 = e v c 2 dt 2 - e x dr 2 - r 2 (d9 2 + sin 2 0d(/) 2 ). (1) 

The pressure and the density satisfy the general relativistic TOV equations of hydrostatic equilibrium |26|,^7| : 

+ (3) 



c 2 r 



dPy = c ( Pu + P v lc 2 ){m + A-Kr 3 P v /c 2 ) 
dr r(r — 2Gm/c 2 ) 

dm 9 , . 

— = 4arr*p u (r), (5) 

where m(r) is the mass enclosed within a radius r. The relevant boundary conditions are m(0) = 0, P W {R) = 0, and 
p v (R) = 0, as the pressure and the density vanish at the radius R of the star. Outside the star, the functions v and 
A are determined by the usual Schwarzschild solution 



e v = e"\ e A 



(l-2GM/c 2 r)-\ (6) 



rR 

M = I Anp v {r)r 2 dr. (7) 
o 
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We now introduce the equation of state of a degenerate relativistic Fermi gas |E8| which may be parameterized as 



P» = K 



X(l + X 2 ) 1 ' 2 [ -X 2 - 1 ) + log [X + (1 + X 2 ) 1 / 2 



(8) 



K 

72 



X(l + X 2 ) 1/2 {2X 2 + 1) - log [X+ (1 + X 2 ) 1 



1/2 



(9) 



8KX 3 



3to„c 2 

Here, n v denotes the neutrino-number density, and K and X are given by 



K = 



X = 



Vv 



(10) 



(11) 



where p„ stands for the local Fermi momentum of the neutrinos of mass m„ , and g u is the spin degeneracy factor of 
neutrinos and antineutrinos, i.e., g v = 2 for Majorana and g v = 4 for Dirac neutrinos and antineutrinos. Using (|^) 
and (H), and introducing dimensionless variables x — r j a v and [i — mjb v with the scales 



av = 2,nH( Mil] 2 Lpi = 2.88233 x 10 10 g^/ 2 
V 9v \ m v J 



17.2 keV 

m v c 2 



km, 



(12) 



= 2J JL ( MziY Mpl = L95197 x 10^ MQg -i/2 ^.2 kcX 



g v \ m v 



m v c* 



(13) 



where M P1 = (hc/G) 1 / 2 and L P1 = {HG/c 3 ) 1 / 2 denote Planck's mass and length, respectively, the TOV equations 
and (|^) can be written as 



dx 



1 + X< 



X(x 2 - 2fix) 



X(l + X 2 ) 1 ' 2 [ -X 2 - 1 ) + log [X + (1 + X 2 ) 1 / 2 



(14) 



— - x 2 
dx 



X{1 + X 2 ) 1/2 {2X 2 + 1) - log fx + (1 + X 2 ) 1/2 



(15) 



subject to the boundary conditions X(0) = Xq and /i(0) = 0. In addition to (14) and (|l5|), there is also an equation 
governing the number of neutrinos n within a radius r — a v x: 



dn 
dx 



= x 2 X 3 {l~2n/x)- 1 ' 2 



where h = h/Nq is the rescaled neutrino-number density subject to the boundary condition n(0) ~ 0, with 



Nn = 



= 3.3765 x 10 
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17.2 keV 

m„c 2 



,-l/2 



(16) 



(17) 



Equations may be solved numerically. Picking up a value Xq for the Fermi momentum at the center 

(in units of mZc), one obtains the total mass of the star M, the radius R. and the total number of particles N, by 
integrating outward until X vanishes. The results are summarized in Figs. |l| and In Fig. [j]the total mass is plotted 
against the radius of the neutrino star. The curve has a maximum, namely, the Oppenheimer-Volkov (OV) limit p7[ , 
at /^ov — 0.15329, which corresponds to a neutrino star mass of 



A I, 



ov 



0.15329 b v = 0.54195 M|im- 2 5i ; 1/2 = 2.9924 x 1O 9 M ( 17 ' 2 k ^ V ^j g -i/2_ 



(18) 



Owing to their large mass, neutrino stars could serve as candidates for supermassive compact dark objects observed 
in the mass range 
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2.5 x 10 6 M© ~ M ~ 3 x 10 9 M (19) 

at the centers of a number of galaxies. Assuming that the most massive and violent objects are neutrino stars at the 
OV limit with Mov = (3-2 ± 0.9) x 10 9 M Q , such as the supermassive compact dark object at the center of M87 M, 
the neutrino mass required for this scenario is 

12.4keV/c 2 < m v < 16.5keV/c 2 for g v = 2, 

10.4keV/c 2 < m„< 13.9keV/c 2 for g v = 4. (20) 

The radius of such a neutrino star is Rov = 4.4466 Rq V , where R ov — 2GMov/c 2 is the Schwarzschild radius of 
the mass Mov- Thus, at a distance of a few Schwarzschild radii away from the supermassive object, there is little 
difference between a neutrino star at the OV limit and a black hole, in particular since the last stable orbit around 
a black hole already has a radius of 3 R v- ^ neutrino star of mass Mov = 3 x 10 9 M Q would have a radius 
Rov = 3.9396 x 10 10 km, or 1.52 light-days. 

Of course, neutrino stars that are well below the OV limit will have a size much larger than black holes of the same 
mass, although they will still be dark and much more compact than any known baryonic object of the same mass. 
As the gravitational potential of such an extended neutrino star is much shallower, significantly less energy will be 
dissipated through accreting matter than in the case of a black hole of the same mass. In fact, there is compact dark 
matter at the center of our galaxy with (2.45 ± 0.40) x 10 6 M concentrated within a radius smaller than 0.0254 pc 
or 30.3 light-days PJlCj], determined from the motion of stars in the vicinity of Sgr A*. Interpreting this supermassive 
compact dark object in terms of a degenerate neutrino star of 2.5 x 10 6 M Q , the upper limit for the size of the object 
provides us with a lower limit for the neutrino mass, i.e., 

m„> 14.3kcV/c 2 for g v = 2, 

m„> 12.0kcV/c 2 for 9lJ = 4. (21) 

In this context, it is important to note that if Sgr A* is a matter-accreting neutrino star p^-p^|, one can, in a 
natural way, explain the so-called "blackness problem" of Sgr A* , i.e., the fact that Sgr A* does not seem to emit 
detectable x-rays of a few tens of keV, which would be emitted by baryonic matter falling towards a black hole. As 
this unmistakable black-hole signature is missing, the concept of a "black hole on starvation" has been created in 
order to save the black-hole idea. However, the neutrino-star model also fits the enigmatic radio-emission spectrum 
of Sgr A* much better than the "black hole on starvation" model pjfl . 

The total mass of the neutrino star M is plotted against the total number of particles N in Fig. |[ For masses 
much smaller than the OV limit, the relation between M and N is unique. However, as M approaches the OV limit, 
M becomes a multivalued function of N . The part of the curve on the left side of the maximum in Fig. |l|, which 
corresponds to the upper part of the curve in Fig. represents unstable configurations p5|p9[ | for which the relative 
mass defect 

Nm u -M , , 

A = TT 22 

eventually becomes negative, as seen in Fig. [|. Thus, for A < 0, the system can gain energy by disintegrating. The 
maximal relative mass defect, or the strongest binding, is obtained at the OV limit with Aov = 3.5807 x 10 -2 . 
For completeness, we note that in the Newtonian limit Xq <C 1, the TOV equations ( |l4| ) and (JlJ) reduce to 

dx x 2 X ' ^ ^ 



ax 6 

which, using the substitution O = X 2 and £ = Ax/ can be cast into the nonlinear Lane-Emden differential equation 
with the polytropic index 3/2 |3(]] 

Owing to the scaling property of the Lane-Emden equation, the mass and radius of a nonrelativistic neutrino star 
scale as Eql 
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In the limit Xq <C 1, (M and (0) yield the equation of state of a nonrelativistic degenerate Fermi gas, i.e., 

/ fi \ 2/3 -4/3*2 

as expected. 

For large central densities Xq ^> 1, /i oscillates around /i,^ = 0.09196, which corresponds to a neutrino star mass 
Moo — 1.795 x 1O 9 M05j / for a neutrino mass m v = 17.2 keV/c 2 . For a gas of neutralinos of a mass precisely equal to 
the neutron mass m n — 0.93955 GeV/c 2 and a degeneracy factor g n = 2, the infinite density limit is M M = 0.4164 M@, 
whereas the OV limit is Mov — 0.7091 M Q and i?ov = 9.1816 km [^9|. Thus, owing to their compactness, neutralino 
stars could mimic the properties of "machos" which have been detected in the dark halo of our galaxy, and which are 
usually assumed to be baryonic brown dwarfs. For large X, the solutions of the TOV equations (|l4| ) and ( plf ) tend to 

3 , „ /3^ 1/4 



The pressure and the density thus become 



and X= - x- 1 / 2 . (28) 



c 4 1 3c 2 1 

Pv = 3 and P" = kr~ "a ' ( 29 ) 
567T r z 5D7T r z 



yielding the equation of state of radiation 



as expected. 



Pu = \<?Pv, (30) 



III. DEGENERATE NEUTRINO AND NEUTRALINO MATTER 

We now turn to the discussion of an astrophysical system consisting of degenerate heavy-neutrino and neutralino 
matter that is gravitationally coupled. As each component satisfies the equation of hydrostatic equilibrium separately, 
i.e., Eq. (pi) and 



dP n 1 2 dv 
— =--(p„c +P n )-, 



the total pressure P = P n + P v and the total mass density p = p n + p v will also obey the same equation 

dP _ 1 .j _|_ p% ^ 
dr 2 dr 



(31) 



(32) 



In addition to the equation of state for neutrino matter, (@) and (Q), we now have the equation of state for neutralino 
matter: 



Pn = K 



9n ( m 



9u \m u 



Y(i + y 2 ) 1 ' 2 ( |r 2 - i J + log (y + (i + y 2 ) 1 / 2 ) 



(33) 



K g n ( m n 

Pn = —— 



Y(l + Y 2 ) 1/2 (2Y 2 + 1) - log [Y+(l + Y 2 ) 1 



1/2 



(34) 



where g n is the spin-degeneracy factor for neutralinos and antineutralinos, and Y is the local Fermi momentum of 
neutralino matter (in units of m„c). Inserting (p3) and (p4j) into (|3l|), after integration we arrive at 



Y =[{\ + Y 2 Y^-^ -I] 1 ' 2 , 



(35) 
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with Yq = Y(0). Using (||), (||), and the equation of hydrostatic equilibrium (0), a similar relation for the Fermi 
momentum of neutrinos (in units of m„c) is obtained: 



X =[(1 + X 2 )c v ^-^ -l] 1 / 2 . 



Combining (35) and (pq), the two local Fermi momenta are related by 

(X 2 + l)Y 2 +X 2 -Y 2 



X 1 



1 + YZ 



The condition X 2 > restricts the range of allowed values of Y to 



Y z > 



YJ-X 2 

i + x 2 



The total pressure and mass density is given by 

P{Y)=P n (Y) + P v {X{Y)) 

and 

p(Y) = p n (Y) + Pv {X{Y)), 



(36) 
(37) 

(38) 

(39) 
(40) 



respectively. 

We now formulate the coupled differential equations describing a gravitationally interacting system of degenerate 
heavy-neutrino and neutralino matter. We first keep the mass of the neutrino halo constant while varying the mass 
of the neutralino star. Introducing the dimcnsionlcss variables x = r/a n and fi = m/b n with the scales 



ml V g„cG 



and br, 



'7T^ 3 C 3 



m l V 9nG 3 ' 



(41) 



the relevant TOV equations can be written in the form 



dY 
dx 



14 


-Y 2 


Y{x 2 


— 2/xx) 




\ 4 


( m u 






) 9n 



Y(l + Y 2 ) 1 ' 2 f^Y 2 - l) + log (Y + (1 + r 2 ) 1 / 2 ) 



[i + x 



X(l + X 2 ) 1 ' 2 [ -X 2 - 1 ) + log [X + (1 + X 2 ) 1 / 2 



(42) 



dp 

dx 



2 {Y(l + Y 2 ) l ' 2 (2Y 2 + 1) + log (Y + (1 + Y 2 ) 1 ' 2 ) 



flu 



X(l + X 2 ) 1 / 2 (2X 2 + l) + log(x + (l + X 2 ) 1 / 2 ) } , (43) 



where X is related to Y through (|37]). If the condition ( p8[ ) is not fulfilled, i.e., the neutrino pressure and density have 
already vanished, the system is solved with the Y terms describing the neutralinos only. 

In order to solve Eqs. ( f42| ) and ( p!3| ) numerically, we fix the Fermi momentum of neutrinos (in units of m v c) at 
the center and vary the central values of the corresponding quantity Yq for neutralinos. The total mass (including 
neutrinos and neutralinos) enclosed within the radius R n of the neutralino star is shown in Fig. ^. Here, the neutrino 
mass and the degeneracy factor are taken to be m v — 17.2 keV/c 2 and g v = 2, respectively, while for the neutralino 
mass we have chosen m n = 939.55 MeV/c 2 and g n = 2, with the scales a n — 6.8304 km and b n = 4.6257 Mq. 
For small neutralino-star masses, the total mass enclosed in R n scales as R^, corresponding to a constant density 
governed by the gravitational potential of the surrounding supermassive neutrino halo. However, as the radius of the 
neutralino star approaches that of a "free" neutralino star, the gravitational potential of the neutralino star becomes 
dominant and the mass now scales as R~ 3 up to the OV limit. Thus there is always a maximal radius of a neutralino 
star within a neutrino halo of a given mass. Substituting neutralinos by neutrons, we must take care of the fact that 
(i) the neutron interacts strongly in the nuclear medium (simulated, e.g., by an effective mass) and (ii) the neutron 
decays through weak interactions. Thus, stable neutron stars can exist only in the range from 0.2 M© to 2 Mq 
where the binding energy is larger than the Q value for the neutron decay. 
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It is instructive to study the properties of a degenerate gas of neutralinos and neutrinos in the nonrelativistic 
approximation. In the limits 1 C 1 and Y -C 1, (|42| ) and ( |43|) simplify to 



dY 
dx 



~x 2 y' 



(44) 



with the boundary conditions 



= |^^ 3 + — f— ) (y 2 + x 2 -y 2 ) 3 / 2 ], 

ax 3 \mr. 



m(o) = o ; y 2 >y 2 -a 2 ; y(o) = y 



(45) 



(46) 



This system of equations can be rewritten in the form of a Lane-Emden type equation by introducing 0„ = Y 2 , 
0„ = X 2 , and a new radial dimensionless radial variable £ = 4x/v / 3 



l__d_ / 2 dQ n 
edZV d£ 



[®T + -( — ) (e„ + e, -e n0 ) 3 / 2 ] , 

5" \ m n 



(47) 



where O„o and O„o ar e the central values of the neutralino and neutrino densities, respectively. For very small 
neutralino densities, i.e., 7<1 and Yq -C 1, the mass equation ( f45|) can be integrated to give 



(48) 



which confirms the conclusion drawn in the context of Fig. [|. 

We now turn to the case of a neutralino star of constant mass surrounded by a neutrino halo of variable mass. The 
TOV equations written in terms of the functions X and fi may be obtained from fl42| ) and ([l3|) , in which we make 
the replacements X <-> Y, g v <-> g n , and to„ <-> m n . Thus, we find 



dA 



A(x 2 - 

4 



/it + x' 



X(l + X 2 ) 1 / 2 ( |X 2 - 1 ) + log (X + (1 + X 2 ) 1 / 2 



^ ^ y(i + y 2 )V2 f F 2_ 1 ) +log ( r + (1 + r2) i /2 

?7i„ J g v \ \6 



(49) 



d/i 
dx 





+ A 2 ) 




4 








5^ 



1/2 



(50) 



with X and Y subject to the condition 



v2 

X 2 >±^ 



Y 2 



Y 2 



(51) 



If this condition is not satisfied, i.e., the pressure and density of neutralinos have already vanished, ([l9]) and (^0|) are 
solved without the Y terms, i.e., for neutrinos only. Choosing the OV limit as the mass of the neutralino star, i.e., 
Mq V = 0.7091 M© for m n = 0.93955 GeV/c 2 and g n = 2, and varying the central Fermi momentum Xo, one can 
find the total mass (including neutralinos and neutrinos) as a function of the radius R v of the neutrino halo. This 
scenario is reflected in Fig. ^ where the length and mass scales are a v = 2.0381 x 10 10 km and b v = 1.3803 x 10 10 M Q , 
respectively. Here the neutrino mass and the degeneracy factor have been chosen as m„ = 17.2 keV/c 2 and g v = 2, 
respectively. At the turning point A, the total mass enclosed within the radius Ra — Roy ~ 9-1816 km of the neutrino 
halo is Ma = Mq V = 0.7091 M Q . At the turning point B, the total mass enclosed within the radius Rb = 0.9912 pc 
of the neutrino halo is Mb = 3.3453 M@. It is interesting to note that, also in this case, there is a maximal radius 
Rb of the neutrino halo, for a given mass of the neutralino star. 

Replacing the neutralino star by a baryonic star, such as a neutron star, a white dwarf, or an ordinary star, the 
only thing that will change in Fig. |5| is the point A at which the enclosed mass starts deviating from the constant 
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value, which depends, of course, on the mass M n of the central object. Thus for M n ~ Mq, the halo will have a size 
of a few light-years and a mass of a few times that of the central baryonic or nonbaryonic star. In this context, it is 
important to note that if every baryonic star is surrounded by such a neutrino halo, the degeneracy pressure of the 
neutrino halo would prevent stars from approaching each other closer than a distance of a few light-years. In such a 
scenario, a large fraction of galactic dark matter would be nonbaryonic. A further attractive feature of this scenario is 
that a neutrino mass of the order of 14 or 15 keV/c 2 , could, at the same time, set the mass scale of the supermassive 
compact dark objects at the centers of galaxies and the scale of interstellar distances in galaxies. 

To investigate the consequences of this idea in more detail, let us assume that the sun is surrounded by a degenerate 
neutrino halo. In the vicinity of the sun, in the region of the size of the planetary system, the neutrino density is 
governed by the gravitational potential of the sun. In fact, the mass due to neutrinos contained within a radius r is, 
in the vicinity of a baryonic or a nonbaryonic star of mass M n , given in the nonrelativistic approximation p^j2^ | by 

^-134xl0 8 a ( M A 3/2 ( Vr r Y /2 (52) 

where AU = 1.496 x 10 8 km is the astronomical unit. This means that for m„c 2 = 17.2 keV, g v — 2, and M n = Mq, 
the mass of the neutrino (and antineutrino) halo contained within the earth's orbit would be M v = 2.68 x 10 _8 M©. 

From the Pioneer 10 and 11 and the Voyager 1 and 2 ranging data J33| we know that the dark mass contained 
within Jupiter's orbit is M d = (0.12 ±0.027) x 10~ 6 M Q and within Neptun's orbit M d < 3 x 10~ 6 Af Q . Of course, the 
Jupiter data should be taken only as a lower limit, as Jupiter tends to eject almost any matter within its orbit Q. 
Nevertheless, taking the Jupiter data at face value, and interpreting dark matter as degenerate neutrino matter, the 
neutrino mass limits are 

12.6keV/c 2 < m v < 14.2keV/c 2 for g y = 2, 

10.6keV/c 2 <m v < 12.0keV/c 2 for g v = 4. (53) 
For dark matter within Neptun's orbit, the neutrino mass limits are 

m v < 15.6keV/c 2 for g v = 2, 

m v < 13.1 kcV/c 2 for g v = 4. (54) 
In summary, considering (p0f), ([H|), and (EjJ), a neutrino mass-range 



14.3keV/c 2 < m v < 15.6keV/c 2 for g v = 2, 

12.0 keV/c 2 < m v < 13.1 keV/c 2 for g u = 4, (55) 



seems to be consistent with all reliable data. 



IV. CONCLUSIONS 



We have studied degenerate fermion stars, consisting of massive neutrinos or neutralinos, or both. We have shown 
that the existence of such objects may have important astrophysical implications. 

For neutrino masses in the range of several keV, neutrino stars are natural candidates for the supermassive dark 
objects at the centers of galaxies. Assuming that the most massive object, such as the compact dark object at the 
center of M87, is a neutrino star at the OV limit, the neutrino mass required for this scenario should be between 
10 keV/c 2 and 16 keV/c 2 , depending on the degeneracy factor g v . 

Furthermore, interpreting the supermassive dark object at the center of our galaxy as a neutrino star, we obtain 
from the upper limit of the size of this object, a lower bound on the neutrino mass which overlaps with the range 
mentioned above. In addition, our interpretation explains the so-called "blackness problem" of Sgr A* in a natural 
way. 

By studying a two-component system consisting of neutralinos in the GeV mass range and neutrinos in the keV 
mass range, we have found that there is always a maximal mass and radius of a neutralino star within a neutrino halo 
of a given mass. Owing to their compactness, neutralino stars could mimic the properties of "machos" . 

Finally, assuming that ordinary stars are surrounded by degenerate neutrino halos of maximal size, a neutrino mass 
of the order of 15 keV/c 2 sets the scale of interstellar distances to a few light-years. 
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FIG. 1. The total mass M of a neutrino star as a function of its radius R. 



FIG. 2. The total mass M of a neutrino star as a function of its total number of particles N. 



FIG. 3. The relative mass defect A as a function of the radius R of the neutrino star. 

FIG. 4. The total mass (including neutralinos and neutrinos) enclosed within the radius R n of the neutralino star for various 
masses M v of the neutrino halo. 

FIG. 5. The total mass of neutralinos and neutrinos contained within the radius R v of the neutrino halo around a neutralino 
star. 
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